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Magnetic  Properties  of  CoCl2-Intercalated  Graphite 

M.  Elahy,*^  M.  Shayegan,^  K.Y.  Szeto*^  and  G.  Dresselhaus^ 
Massachusetts  Institute  of  Technology 
ibridge,  MA  02139,  USA 


The  temperature  and  magnetic  field  dependence  of  the  magnetic 
susceptibility  and  heat  capacity  of  graphite-CoCl^p are  investigated  to 
understand  the  role  of  dimensionality  in  the  observed  magnetic  phase 
transitions.  The  structure  of  graphite-CoClo1 shows  similarities  between  the 
magnetic  layers  of  Co*2  in  pristine  and  in  the  Intercalated  CoCl^ compound . 
However,  micros tructural  analysis  shows  the  presence  of  small  islands  of 
magnetic  ions.  The  theoretical  model  for  the  2-dimensional  planar  magnets  is 
extended  to  include  small  magnetic  domains.  The  finite  size  effect  analysis 
shows  good  agreement  between  the  experimental  and  theoretical  results.  The 
presence  of  external  magnetic  fields  demonstrates  the  static  scaling 
hypothesis  for  the  magnetic  susceptibility  and  heat  capacity. 
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INTRODUCTION 

Intercalation  of  magnetic  species  Into  graphite  provides  an  Ideal  system 
for  the  study  of  two-dimensional  magnetism  and  for  the  investigation  of  the 
role  of  dimensionality  in  magnetic  and  structural  properties.  In  these 
systems  the  interplanar  interaction  of  magnetic  spins  can  be  systematically 
reduced  by  Inserting  a  controlled  number,  n,  of  diamagnetic  graphite  layers 
between  each  adjacent  pair  of  magnetic  Intercalate  layers,  where  n  is.  the 
stage  index  of  the  graphite  Intercalation  compound  (GIC). 

Recent  susceptibility  measurements  on  FeCl3-GIC,  C0&2-GIC  and  NICI2-GIC 
have  shown  similarities  in  both  the  temperature  and  the  magnetic  field 
dependence  of  their  magnetic  behaviors  [1].  It  has  furthermore  been  shown  [2] 
that  for  these  compounds  the  fit  of  the  susceptibility  above  the  critical 
region  to  a  classical  Ourle-Ueiss  law  is  poor.  A  much  better  fit  is  obtained 
using  the  results  of  the  Kostcrlitz-Thouless  theory  (K-T)  for  the 
two-dimensional  xy  model  [3]. 

We  have  chosen  the  Co&2~G1C  system  for  a  more  critical  investigation  of 
this  phase  transition  because  it  has  the  largest  xy  anisotropy  in  the 
spin-spin  Interaction  among  the  known  magnetic  graphite  Intercalation 
cosipounds . 

The  lattice  structure  of  pristine  C0CI2  is  trigonal  with  space  group  R3m, 
consisting  of  hexagonal  layers  of  Co+^  in  a  triangular  net  separated  from  each 
other  by  two  hexagonal  layers  of  Cl".  Pristine  C0CI2  undergoes  an 
antiferromagnetic  phaae  transition  at  the  Neel  temperature  Tjj  *  24.9  K.  For 
temperatures  below  Tg,  the  spins  lie  in  the  layer  plane  and  align  in 
two-dimensional  ferromagnetic  sheets  which  are  stacked  antiferromagnetically. 
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Experimental  modeling  of  the  neutron  scattering  experiments  [4]  have 
established  values  for  the  intraplanar  and  interplanar  exchange  energies  of 
J  •  28*5  K  and  J'  ■  2.2  K,  with  c-axis  anisotropy  energies  of  D  *  -16.0  K  and 
D*  -  -3*3  K,  respectively.  The  in-plane  anisotropy  field  is  estimated  to 
be  HA  -  1.0  x  10~2  x  from  the  spin  flop  field  Hgp  ■  0.2  T,  and  the 
saturation  field  Hg  •  3.3  T  [5]. 

This  paper  focusses  on  two  important  aspects  of  magnetic  phase 
transitions  in  these  magnetic  61C.  The  first  topic  concerns  the  finite  size 
of  the  magnetic  domains  and  its  effect  on  the  observed  transition  phenomena. 
The  second  topic  concerns  the  role  of  an  external  magnetic  field  on  the 
observed  magnetic  phenomena. 

SAMPLE  PREPARATION  AND  CHARACTERIZATION 

Samples  of  graphlte-CoCl2  vers  prepared  by  the  standard  two-zone  vapor 
transport  method  for  graphite  intercalation  compounds  [6,7].  Prior  to  heat 
treatment,  the  pristine  C0CI2  and  the  host  HOPG  sample  were  separated  at  a 
distance  of  a  few  cm  inside  a  pyrex  tube.  Before  sealing  the  ampoule,  an 
atmosphere  of  CI2  (at  ~  350  torr  pressure)  and  some  AICI3  (~  5Z  of  C0CI2  by 
weight)  were  added  to  the  ampoule  to  act  as  a  complexing  halide  [6].  The 
presence  of  this  AICI3  was  found  to  be  crucial  to  the  preparation  of 
well-staged  samples  for  n  >  2.  A  two-zone  furnace  was  used  with  the  C0CI2 
heated  to  ~  500  K  and  the  HOPG  to  a  slightly  higher  temperature.  A  larger 
temperature  difference  between  the  graphite  and  the  C0CI2  zones  results  in  a 
higher  stage  of  the  Intercalation  compound.  The  intercalated  samples  were 
characterized  for  stage  index  and  stage  fidelity  by  (00!)  x-ray 
dif f r aerograms . 
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To  establish  the  role  of  the  AICI3  complexing  halide,  we  have  done 
extensive  structural  characterizations  of  our  samples.  The  results  of  these 
characterizations  are  presented  here  and  are  applied  to  the  interpretation  of 
the  results  of  our  magnetic  measurements.  The  microstructure  was  examined 
using  a  scanning  transmission  electron  microscope  (STEM)  and  the  presence  of 
▲ICI3  was  investigated  by  x-ray  fluorescence  techniques. 

This  analysis  shows  that  AICI3  and  C0&2  do  not  mix  in  the  Intercalated 
compound  and  each  species  forms  its  own  domain.  In  Fig.  1(a)  a  typical  (STEM) 
micrograph  of  the  boundary  of  C0&2  and  AICI3  is  shown.  The  darker  circular 
areas  are  AICI3  re8lon*  and  the  bright  parts  are  C0&2  domains.  Within  the 
resolution  of  the  x-ray  fluorescence  instrument,  no  A1  was  observed  in  the 
C0&2  domain  and  no  Co  in  the  AICI3  domain.  Thus  our  system  is  modeled  by 
small  crystallites  of  staged,  ordered  CoC^-GIC.  The  measurement  of  the 
volume  of  intercalated  C0CI2  to  the  volume  of  Intercalated  AICI3  in  the  GIC 
samples  is  included  in  the  analysis  of  our  magnetic  measurements.  These 
stoichiometric  determinations  are  done  by  chemical  analysis  of  the  bulk 
samples. 

Figure  1(b)  shows  an  electron  diffractogram  for  the  same  sample  as  in 
Fig.  1(a).  The  diffractogram  is  Interpreted  as  a  superposition  of  pristine 
graphite  spots,  pristine  C0CI2  spots  and  rings  whose  radii  match  the  lattice 
constants  of  AICI3  (see  Fig.  1(c)).  Thus  the  electron  diffraction  results 
further  confirm  the  separation  of  AICI3  and  C0CI2  into  distinct  domains.  The 
results  further  show  an  angular  correlation  between  the  C0CI2  layer  and  the 
adjacent  graphite,  layers.  The  electron  diffraction  pattern  implies  a 
completely  random  angular  distribution  of  the  AICI3  planes,  and  Indicates  that 
the  structure  of  C0CI2  In  C0CI2-GIC  is  the  same  as  in  its  pristine  form.  This 
experimental  result  indicates  that  the  known  in-plane  magnetic  parameters  of 
pristine  C0CI2  and  C0CI2-GIC  are  naarly  identical. 
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We  also  examined  Che  structure  of  Che  C0CI2  domains  using  che  dark  field 
Image  Cechnlque  and  found  chac  Chese  domains  consisc  of  small  (100  A  x  500  A) 
cigar  shape  Islands  which  are  uniaxi ally-oriented  and  concaining  -9000  spins. 

•  Figure  1(d)  shows  a  typical  dark  field  image  of  C0CI2  domains.  The  cigar 
shape  of  Che  domains  is  aCCribuCed  Co  a  minimization  of  the  energy  by  aligning 
Che  maximum  number  of  incerealace  atoms  over  che  hollows  of  che  graphite 
honeycomb  structure.  The  uniaxial  orientation  of  che  C0CI2  domains  allows  for 
a  high  in-plane  and  c-axis  packing  density.  Although  Che  C0CI2  remains 
incommensurate  with  Che  graphite  substrate,  this  oriented  needle-shaped  domain 
arrangement  creates  an  angular  locking  between  the  two  lattices.  This  locking 
is  also  clear  in  the  diffraction  picture  (Fig.  1(b)).  The  presence  of  these 
intercalant  Islands  requires  an  extension  of  present  2-dimensional  theories  to 
include  finite  size  effects,  as  discussed  below. 

In  order  to  quantify  our  magnetic  measurements,  the  chemical  formulae 
CgaX,  where  X  stands  for  (CoC^i.y  (AlC^y  were  determined  from  chemical 
analysis  [8] .  Using  our  (00i)  x-ray  stage  determination,  the  parameters  £  and 
y  were  found  to  be  7.07  and  0.68  respectively  for  the  stage  4  compound.  All 
compounds  measured  for  stages  n  ■  2,  3,  4,  and  5  were  found  to  yield  £  and  y 
values  in  the  range  7  <  £  <  8  and  0.63  <  y  <  0.73)  [8] . 

RESULTS  AMD  DISCUSSION 

In  Fig.  2a  we  plot  our  in-plane  susceptibility  results  for  stage  4 
C0CI2-GIC  in  the  absence  of  an  externally  applied  magnetic  field.  The  shape 
of  the  temperature  dependence  has  been  shown  to  fit  a  two-dimensional  model  of 
the  Kosterlitx-Thouless  form  [3].  The  applicability  of  such  a  2-D  model  to 
our  system  can  be  understood  from  the  following  considerations. 
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On  Che  basis  of  our  structural  studies  of  graphite-CoCl2,  we  can  assuse 
the  following  magnetic  Hamiltonian  for  the  system  in  the  absence  of  an 
external  magnetic  field: 


H  -  J  2  H  Sj  +  JX  Z  Sjz  +  J*  Z  Si  Sk 

ij  ij  ij 


where  the  values  of  the  intraplane  coupling  constant  J,  and  its  associated 
anisotropy  JA  are  taken  to  be  essentially  the  same  as  in  the  pristine  C0CI2 
compound,  i.e.  J  -  28.5  K  and  JA  ■  -16.0  K,  respectively,  the  introduction  of 
several  diamagnetic  graphite  layers  is  expected  to  further  decrease  the 
already  small  value  of  J'  in  the  pristine  compound.  In  fact,  our  magnetic 
field  dependent  behavior  of  the  susceptibility  indicates  that  J'/J  in  a 
stage  2  compound  has  an  upper  bound  of  2  x  10~*.  The  dipole-dipole 
interaction  for  two  spins  separated  by  the  ad3acent  magnetic  layer  distance  it 
stage  2  gives  a  lower  limit  for  this  ratio  o'f  ~  6  x  10“5  K. 

In  the  high  temperature  phase  characterized  by  unbound  vortices,  the  K-T 
correlation  length,  5,  increases  with  decreasing  temperature.  When  £  grows  ti 
~  Sq  (J/JA)  ■  1.3  a0,  where  a0  is  the  lattice  constant,  the  energy  required  t 
rotate  a  spin  in  the  direction  of  the  hard  axis  (c-axis)  becomes  of  the  same 
magnitude  as  J.  thus,  below  28.5  K,  where  the  thermal  energy  of  the  spins 
becomes  less  than  the  spin-spin  coupling,  we  expect  the  spins  to  favor  an  xy 
system.  On  the  other  hand,  at  much  lower  temperatures, 

e/e0  ~  (J/J')*^  ~  100,  and  the  effect  of  the  adjacent  magnetic  plane  is  felt 
If  we  assume  that,  on  the  average,  each  island  contains  less  than  104  spins, 
then  the  system  should  satisfy  the  2-D  criterion  throughout  the  critical 
region.  Ma  thus  conclude  that,  in  the  temperature  range  of  interest,  the  2-1 
xy  model  is  appropriate  to  the  graphlte-CoCl2  system. 
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The  Kosterlitz-Thouless  theory  for  the  two-dimensional  xy  model  has  been 
extended  by  Jose,  Kadanoff,  Kirkpatrick  and  Nelson  ( JKKN)  to  include  the 
effect  of  the  six-fold  symmetry  breaking  anisotropy  field  [10].  They  found 
that  there  exists  a  critical  region  in  temperature  with  infinite  initial 
susceptibility  on  both  sides  of  which  x  diverges  as  exp(b|t|-1/2).  The  upper 
temperature  divergence  of  the  susceptibility  corresponds  to  a  K-T  type  vortex 
binding-unbinding  phase  transition  while  the  lower  temperature  divergence 
corresponds  to  a  transition  from  the  bound  vortex  state  to  a  state  with 
conventional  long-range  order. 

The  analysis  of  JKXN  Is  valid  only  for  t  <  0.1  where  t  -  T/Tc  -  1.  But 
we  have  extended  the  JXXN  model  to  temperatures  further  above  the  upper  phase 
transition  temperature  Tc}  using  the  Nelson-Rudnick  extension  formalism  and 
the  high  temperature  2-D  Coulomb  gas  Green  function.  We  found  that  the  same 
functional  form  is  applicable  to  a  much  larger  temperature  range,  with  the 
parameter  b  being  weakly  temperature  dependent  [9] . 

Thus  the  use  of  the  exponential  function  to  fit  our  experimental  data 
over  a  wide  temperature  range  as  reported  earlier  can  be  theoretically 
justified.  The  same  experimental  situation  has  been  reported  for  K2CUF4  [11]. 
We  have  also  found,  using  numerical  computation,  that  the  rounding  of  x  the 
susceptibility  peak  is  consistent  with  a  finite  size  analysis,  using  values 
for  the  island  size  found  experimentally. 

As  reported  previously  [12]  the  magnetic  heat  capacity  measurements  are 
consistent  with  the  Monte  Carlo  calculations  based  on  the  xy  model  [13]. 

Both  experimental  and  theoretical  results  exhibit  a  broad  peak  slightly  above 
the  temperature  of  the  phase  transition  from  unbound  to  bound  vortices. 
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The  application  of  an  external  magnetic  field  (in-plane)  has  a  dramatic 
effect  on  the  magnetic  properties.  In  the  case  of  the  susceptibility,  small 
magnetic  fields  have  significant  effects,  as  shown  in  Fig.  2b  where  we  plot,  /. 
vs.  temperature  for  several  values  of  the  externally  applied  magnetic  field 
H.  As  H  increases,  the  magnitude  of  the  susceptibility  at  the  peak, 
decreases,  and  the  position  of  the  peak  in  temperature,  Tmjv(X),  shifts  to 
higher  temperatures. 

In  Fig.  3  we  plot  the  field  dependence  of  L,  and  on  a  log-log  scale. 

The  data  fit  a  power  law.  of  the  form  Xux(H)  ~  HT^  where  \  -  0.80.  Similar 
magnetic  field  and  temperature  dependences  have  been  reported  by  Suzuki  and 
Ikeda  in  K2CUF4  [14],  who  have  found  agreement  between  their  results  and  the 
static  scaling  hypothesis  of  Kadanof f  [15]  •  The  spin-spin  interaction  in 
K2CUF4  is  primarily  of  the  Heisenberg  type  with  only  1%  xy  anisotropy  and  the 
in-plane  symmetry-breaking  field  is  negligible.  Nevertheless,  K2CuF 4  and 
C0CI2-GIC  exhibit  essentially  identical  behavior  in  the  pertinent  T  and  H 
ranges.  The  value  of  \  ■  0.82  is  obtained  for  K2CUF4  [14],  in  excellent 
agreement  with  our  values  for  \  in  CoC12~GIC. 

Both  the  X(T)  and  0^(1)  measurements  on  CoCl2-0IC  show  that  T^aax 
shifts  to  higher  temperatures  with  Increasing  field  as  T^tCH)  -  T^aX(0)  -  h'- 
(see  Fig.  3).  From  this  plot  we  estimate  the  exponent  \i  from  the  magnetic 
field  dependence  of  the  maximum  in  X  to  be  ■  0.37  t  0.05  and  from  to  be 
(t  ■  0.9  ±  0.3.  The  discrepancy  between  the  value  of  \i  obtained  from  the  field 
dependence  of  the  susceptibility  and  of  the  heat  capacity  in  the  present  work 
is  not  fully  understood. 
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Our  results  Indicate  that  in  graphite-CoCl2»  when  the  amplitude  of  the 
external  field  exceeds  the  in-plane  6-fold  symmetry  breaking  fi»ld  of  100  0e, 
the  behavior  of  %  is  also  in  good  agreement  with  the  static  scaling 
hypothesis.  Furthermore,  the  critical  exponents  in  this  magnetic  field  regioi 
are  independent  of  stage  index.  This  suggests  that  the  interaction  of  a  spin 
with  the  external  field  in  this  region  dominates  over  the  interplanar 
spin-spin  interaction.  This  in  turn  suggests  a  maximum  value  of  J'  for  the 
lowest  stage  (n  ■  2)  which  is  J'nax  *  7  x  10-3  K. 


CONCLUSIONS 

The  experimental  CoCl2~GIC  system  provides  an  important  experimental 
system  for  study  of  the  2-dimensional  xy  model.  The  static  scaling  hypothesi 
as  well  as  finite  size  effects  are  introduced  to  interpret  the  measurements. 
This  model  yields  excellent  agreement  with  the  experimental  data. 
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FIGURE  CAPTIONS 

Fig.  l.(a)  STEM  bright  field  image  -The  dark  regions  are  AICI3  domains  and  the 
bright  regions  surrounding  them  are  C0CI2  domains.  X-ray 
fluorescence  shows  no  C0CI2  in  the  AICI3  regions  and  no  AICI3  In  the 
C0CI2  regions. 

(b)  Electron  diffractogram  of  a  graphite-CoCl2  sample.  To  within  our 
resolution,  the  pattern  shows  a  super-position  of  HOPG  and  pristine 
C0CI2  diffraction  spots.  Although  C0CI2  is  incommensurate  with  the 
graphite  substrate,  a  high  degree  of  angular  locking  is  observed 
between  them,  and  hence  the  spot  pattern.  No  such  correlation  is 
apparent  between  the  graphite  and  AICI3  diffraction  peaks  which  form 
circular  rings  with  radii  compatible  with  the  inplane  lattice 
constant  of  pristine  AICI3. 

(c)  A  schematic  representation  of  the  diffraction  pattern  of  Fig.  lb. 

(d)  TEM  dark  field  image  of  C0CI2  diffracted  light  in  graphite-CoCl2> 
The  C0CI2  regions  consist  of  needle  shaped  Islands  of  100A  x  500A 
dimensions,  uniaxlally  oriented  to  allow  for  maximum  packing  and  a 
minimum  in  free  energy. 

Fig.  2.  Temperature  dependence  of  the  in-plane  susceptibility  of  stage  4 

CoCl2~GXC  for  (a)  H  ■  0  and  (b)  for  several  external  magnetic  field 
values. 

Fig.  3.  Field  dependence  of  the  susceptibility  peak  amplitude  X^x  (right 

scale),  and  of  the  shift  in  the  susceptibility  peak  position  T^ffy(x) 
and  the  heat  capacity  peak  position  (right  scale). 
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